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Cavity quantum electrodynamics (CQED) [1] plays an elegant role of studying strong coupling between light 
and matter. However, a non-mechanical, direct and dynamical control of the used mirrors is still unavailable. 
Here we theoretically investigate a novel type of dynamically controllable cavity composed of two atomic mir¬ 
rors. Based on the electromagnetically induced transparency (EIT) [2, 3], the reflectance of atomic mirror is 
highly controllable through its dispersive properties by varying the intensity of applied coupling fields or the 
optical depth of atomic media. To demonstrate the uniqueness of the present cavity, we further show the pos¬ 
sibility of manipulating vacuum-induced diffraction of a binary Bose-Einstein condensate (BEC) [4, 5] when 
loading it into a dispersive cavity and experiencing superradiant scatterings [6]. Our results may provide a novel 
all-optical element for atom optics [7] and shine new light on controlling light-matter interaction. 


For decades, CQED has served as an elegant model for 
demonstrating atom-light interaction and fundamental princi¬ 
ples of quantum mechanics [1]. In conjunction with the state- 
of-art experiments using ultracold atoms, CQED has been able 
to facilitate the realization of quantum metrology [8], quan¬ 
tum teleportation [9] and quantum memory storage [10], etc., 
experiments paving the way for the fulfillment of quantum 
technology. By its very nature, cavity is an indispensable in¬ 
gredient to CQED, which reflects photons and has them in¬ 
teract with atoms in a confined space. However, conventional 
CQED systems consisting of reflecting mirrors limiting them¬ 
selves from being networked at large scale. Here, we propose 
a novel cavity system where mechanical mirrors are replaced 
by atomic ensembles whose dispersion can be controlled by, 
e.g., EIT, and so we term it as dispersive cavity which of¬ 
fers unique controllability. Moreover, the separation between 
the atomic mirrors and their size can be much smaller than 
that of mechanical ones. Therefore, the cavity cooperativ- 
ity is enhanced due to smaller cavity volume and it becomes 
more feasible to reach strong coupling regime and scalable for 
more demanding systems. Both advantages render fabricating 
a more functional matter-wave interferometry or atomic cir¬ 
cuit possible. The dispersive cavity can also be integrated on 
atom chips [11] or used to study the cavity-mediated spin- 
orbit coupling in the ultracold atoms [12], bringing both the 
controllability and the tiny size into full play. 

As depicted in Fig. 1 a, our dispersive cavity is com¬ 
posed of two remote collections of three-level-A type atoms. 
Based on EIT [13, 14], as illustrated in Fig. 1 b and d, each 
atomic medium constitutes a mirror when interacting with two 
counter-propagating coupling fields with Rabi frequencies Q* 
driving the |2) |3) transition. Under the action of four-wave 

mixing [14], an incident probe field with Rabi frequency £l + 
driving the transition |3) <-> |1) will be reflected and generate 
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a backward probe field with Rabi frequency Q~ [13], and vise 
versa. For showing the controllability and novelty of the dis¬ 
persive cavity, we consider an excited quasi one-dimensional 
binary BEC which transversely flies through it and emits pho¬ 
tons [5]. As illustrated in Fig. 1 c, the binary BEC consists 
of atoms condensed in the ground state | g) and excited state 
\e). The photonic polarization and frequency of transition 
| g) —> \e) match that of |1) —> |3) of two atomic mirrors such 
that emitted photons with wavelength A from BEC can be re¬ 
flected by two side EIT media. The vacuum fluctuations [15] 
and the geometric shape of the BEC will speed up the decay 
and result in the emission of superradiance (SR) along BEC’s 
long axis [6, 15, 16]. Therefore when the long axis of the 
loaded BEC is parallel to the cavity axis, the SR will circulate 
within a dispersive cavity and become subsequently releasable 
by switching off coupling fields. The coexistence of BEC and 
superradiant photons within a cavity enables the multiple scat¬ 
tering between them. This renders the generation of matter 
waves with high order harmonic of photon wave number pos¬ 
sible, namely, the diffraction of matter waves [7]. Because 
the effect is triggered by vacuum fluctuations, we term it as 
vacuum induced diffraction of matter waves. 

We first investigate the controllability of the novel cavity 
and calculate the reflectance of an EIT mirror by numerically 
solving the optical-Bloch equations [14, 17] 

dtp = ^ [#,p] +pdec , (1) 

(- 3 , ± = itjEiTp^i , ( 2 ) 

where p is the density matrix for the state vector 2^=1 ^'10 
of the three-level atom in the Fig. 1 a, b and d. Also, pd ec 
describes the spontaneous decay of the excited state |3) char¬ 
acterized by rate T. H is the Hamiltonian for the interaction 
between atomic mirrors and counter-propagating fields (see 
Methods). Moreover, t/eit = Td opt /2L with d opt and L being 
the optical depth and the length of the medium, respectively. 
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FIG. 1. Schematic plot of a dispersive cavity, a, a binary BEC sandwiched by two three-level EIT media, where the former is initially in the 
state \e) and the latters are in state |1). The red arrows represent two counter-propagating coupling fields driving the transition |2) —> |3) as 
shown in b and d. The quantum fluctuations of binary BEC spontaneously generate the forward and backward SR fields Q* (yellow) as shown 
in c, which connect the EIT transition from |1) —» |3). The composition of three-level atoms and coupling fields acts as a dispersive cavity, 
which reflects resonant yellow photons. The controllable reflectance of a dispersive mirror and the enhancement of cavity finesse are shown 
in e and f, respectively, g, the coherently trapped superradiant field [|f2 + (t,z)| 2 + |Q _ (t,z)| 2 ]/T 2 is shown where EIT mirrors are located in the 
regions |z| >0.2 mm. 


The resultant dispersive cavity has only one reflecting band 
whose width is associated to the EIT transparency window 
A qjeit = £l 2 c / (r ^d opt ^ which is highly tunable by either vary¬ 
ing the optical depth or the coupling fields. As shown in Fig. 1 
e, at various detuning and a fixed optical depth d opt = 500, the 
reflectance approaches unity when the fields D* are resonant 
and the width can be controlled by varying the intensity of 
the coupling fields. To gain more insight into the dispersive 
cavity, we invoke the effective cavity finesse /( A p ) = ) 

where R(A p ) is the reflectance as a function of the detuning 
A p of fields D* [18]. The cavity finesse is plotted in Fig. 1 f, 
which demonstrates that the cavity finesse can be controlled 
by tuning the EIT transparency window [2, 3]. 

We are ready to show the cavity effect on the BEC dynamics 
by numerically integrating equations (1) and (2) together with 
the coupled Gross-Pitaevskii equation which reads [19, 20] 
(see methods) 


ihd t 


if/ g 

4>,e 


H sp + H s im H atom -u g ht \ / \j/ g \ . 

H :,om-u gh , Hsp+Kt-'i /Ur (j 


where i// g , e are the wave functions of the binary BEC, H sp the 
single-particle Hamiltonian, H 8 ^ the atom-atom interaction, 
H atom-light the atom-light interaction (see Methods), and T is 
the spontaneous decay rate of the excited level \e) introduced 
phenomenologically. To start up the superradiant process of 
the binary BEC in a cavity, the vacuum fluctuation is intro¬ 
duced by truncated Wigner method [21] (see methods). The 


dynamics of the generated SR fields is described by 

(± Q* = irjBEccrfg^, (4) 

where cr^ is the coherence of the BEC associated with the 
plane-wave factor e ±lk P z (see Methods for details) and rj BEC = 
3 TNbecA 2 /4nA with A being the transverse cross section of 
the condensate. In the following, we consider the EIT mirror 
of optical depth d opt = 500 and coupling field of £l c = 2r. A 
quasi one-dimentional BEC of 7 Li at 2 Si/ 2 (|g» state can be 
prepared in a single beam dipole trap. An additional TEMoi- 
type laser mode beam is applied to cross the atomic ensemble 

[22] . The spatial intensity distribution of the cross beam al¬ 
lows to trap low density thermal clouds aside as the EIT mir¬ 
rors while leaving the central BEC undisturbed. The 7 Li BEC 
can be coherently transferred to the excited 3 P 3/2 (\e)) state by 
a 7r pulse. The narrow linewidth of 3 P 3/2 state, T ~ 760 kHz 

[23] offers a reasonable time scale for detection. Moreover, 
d opt = 1000 has already been experimentally achieved in cold 
atom systems [24, 25]. As an example, with Nbec = 3 x 10 5 
and A = (n x 3 2 ) yum 2 , Fig. 1 g depicts the intensity of the 
confined superradiance \£l + (t,z)\ 2 + \£l~(t,z)\ 2 in a dispersive 
cavity. 

The three steps of vacuum induced diffraction of matter 
waves is shown in Fig. 2 a, namely, from top to bottom, cre¬ 
ating BEC in state \e), loading BEC into a dispersive cavity 
and releasing BEC. The multiple-scattering between circulat¬ 
ing SR fields and the BEC generates higher order harmonics 
of the matter wave. Theses high k p modes can be clearly ob¬ 
served by time-of-flight (TOF) simulation with Nbec - 10 6 
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FIG. 2. Controllable vacuum induced diffraction of matter waves using a dispersive cavity, a, coherent matter waves initially prepared 
in its internal excited state are diffracted by three steps from top to bottom. The condensate is created in a harmonic trap (green surface) and 
flies through a dispersive cavity where odd multiple k p modes of matter waves emerge. When switching off the trap, matter waves of different 
modes split due to the momentum kick from SR photons, namely, atomic diffraction, b, the BEC’s time of flight dynamics. Six density bumps 
carrying momenta of ±k p , +3 k p , and ±5k p can be clearly observed at the end of time, c, the coupling-field-strength-dependent mean particle 
number of different modes of diffracted BEC carrying ±k p (red), +3 k p (green) and ±5k p (blue). All the data points are averaged over 1000 
realizations, and the color filled regions illustrate the corresponding error bars. 




FIG. 3. Dynamical control of vacuum induced diffraction of matter waves using a dispersive cavity, a, when loading BEC into a dispersive 
cavity, atomic diffraction can be dynamically controlled by turning off the coupling field at different instance. The time of flight for that 
coupling field is off at b 0.02 fis and c 0.07 fis after loading BEC. The particle numbers on each mode of b and c are (N ±kp , N ±3kp , N ±5kp ) =(8.6x 
10 5 ,300,0) and (8.2 x 10 5 ,1.6 x 10 4 ,200), respectively. 


as shown in Fig. 2 b. Especially at txoF = 1.5 ms there are 
six notable peaks corresponding to -5 k p , -3 k p , -k p , k p , 3 k p 
and 5k p from left to right. To study the controllability of gen¬ 
erating high-k p modes, we tune the cavity finesse by varying 
the strength of coupling fields and perform 1000 independent 
realizations for statistics. As shown in Fig. 2 c, the population 
of ±k p modes is gradually and coherently transferred to that 


of higher order harmonics while increasing Q c . The particle 
number of ±3k p and ±5 k p modes are significantly generated 
when > O.ir and IT, respectively. This is in contrast to the 
situation without the dispersive cavity where only ±k p modes 
and very small fraction of ±3 k p modes are generated, namely, 
about a factor of 1000 enhancement by the cavity for ±3 k p 
modes. Remarkably, ±5 k p modes only emerge under the ac- 
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tion of a dispersive cavity. 

Figure 3 illustrates the possibility of dynamically manipu¬ 
lating atomic diffraction. In Fig. 3 a we demonstrate a pro¬ 
tocol to suppress the generation of particular harmonic of k p 
matter wave by switching off two coupling fields at some in¬ 
stance, which dynamically disables the dispersive cavity. We 
respectively turn off the coupling fields at 0.02 /is and 0.07 /is 
in Fig. 3 b and Fig. 3 c after BEC is loaded. Our TOF sim¬ 
ulations show that ±3k p and ±5k p modes can be successfully 
suppressed by release SR photons at proper instances before 
corresponding multiple scattering happens. 

In conclusion, the novel flexibility of a dispersive cavity 
may give the possibility of implementing more delicate ma¬ 
nipulations of matter wave and becomes all-optical elements 
for atom optics. Additionally, an all-optical Q-switching su¬ 
perradiant source may result from combining our system with 
triggered superradiance [26]. Moreover, the recent advance 
on x-ray quantum optics [27, 28] also shows the potential of 
achieving atomic mirrors in x-ray domain [29]. Given the 
high, non-mechanical and dynamical controllability of the dis¬ 
persive cavity , the present system provides a unique and novel 
stage of manipulating light-mater interaction. 


METHODS 

In the presence of the counter-propagating coupling and SR 
fields, the Hamiltonian of the EIT medium is H = Ho + Hj 
where Ho = 2 / I j) 01 and Hj is the atom-light interaction 
Hamiltonian of the EIT medium given by form as 

&i = \ |3> <11 

+ (Q+e^cz-^0 + n - e i(-fez-^c*)) | 3 > <2| + H.c] . 


In equation (3), the operators accounting for the quasi ID 
condensate dynamics take the form as H sp = 

^atom-light = § {&* e~ ik ^ + £TW) and H\ nt = gjj \H 1 2 + 
gji 1071 2 which are the single-particle, atom-light interaction, 
and nonlinear interaction Hamiltonians, respectively. The 
quasi ID nonlnear interaction constant takes the form as gji = 


2 7rajiQj r N bec 

(^Z&OSC 


where a osc 



is oscillator length. The oj r and 


aj z are the transverse and longitudinal trapping frequencies, 
respectively, where we use ( co r , co z ) = 2n x (100,2.1 x 10 3 4 5 6 7 8 9 10 11 12 ) Hz 
in the simulation. The optical coherence of the condensate 
is given by cr eg = 0^0*. In the presence of circulating SR 
fields, the macroscopic wave functions and coherence would 
consist of high ±nk p modes, i.e. ij/ g = Yj7=-oo 0 > e 


= 'Zn=-oo'/' ( e ) e i2nkpZ , and CT eg = Zn=-oo ^eg e inkpZ with l//\ 


(n) ^i(2n+l)k p z ^ 

,/n) 

J e,g 

and crf g being the spatially slowly-varying envelope. The co¬ 
herence terms, <xcarrying e ±lk P z account for the generation 
of SR fields. 

The dynamics of the condensate flying through the cavity 
is simulated by the initial state where all the atoms are con¬ 
densed in the excited state forming a Thomas-Fermi density 
profile while the ground state consists of quantum fluctuations 
sampled according to the truncated Wigner approximation of 


= VT7 Nbec 2 "i a£j fe), < a*ai >= 6j,/2, and £ (z) be- 
ing an orthonormal basis. 

Equations (l)-(4) are numerically solved by the method of 
lines where the Bloch and GP equations are propagated by 
the Fourier pseudo spectral method and the adaptive Runge- 
Kutta method of orders 4 and 5 (RK45) for space and time 
integration, respectively, while the SR fields are integrated by 
the semi-Euler method. In the simulation, the total number of 
Fourier modes is 5 x2 12 and the number of mode for sampling 
initial fluctuations is M = 3000. 
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In this supplementary material, we present the equation of motion and the condensate dynamics in more detail. 

In the presence of counter-propagating coupling fields Q* and the superradiant fields Q*, the dynamics of the EIT mirror is 
described by the optical-Bloch equation which reads: 

d t pu = -i (A - A c )pi2 - z'(ffp | 3 + ^cPn )/ 2 + K&*P 32 + ^”*P 3 2)/ 2 > (SI) 


d tP f 3 = - (/A + r/2)p* 3 - if2fp 12 /2 - ift** (pn -P 3 3 ) /2, 


(S2) 


d,p% = - UA C + r/ 2 ) p% - AV*P 21 + inf (p 33 - p 22 ), (S3) 

fpn = rp 33 /2 + im [n + pf + n _ pf ], (S4) 

dtpu = rp 33 /2 + im [off , + n,:p 2 3 |, (S5) 

d,p 33 = -rp 33 - im [n + Pi 3 + n _ pi 3 + pf + p^j, (S 6 ) 


where T is the spontaneous decay rate of the excited state |3) and pij is the element of the density matrix. In equations (Sl)- 
(S 6 ), all fast-oscillation exponential factors associating with center frequencies and wave factors have been eliminated, and only 
slowly-varying profiles are retained. 

When the condensate is loaded in the dispersive cavity, in the presence of the interaction between the circulating SR fields 
and condensate, the ground-state and excited condensate wave functions constitute the superposition of discrete ±( 2 n + 1 )k p and 
+2 nk p plane waves, respectively. In this manner the condensate wave function can be decomposed as : 

oo 

ij/ g = J] ^«y( 2 «+i ( S7) 


and 


oo 

ike = Yj 'k ( ? ) e i2nkpZ ’ 

n=-oo 


where ij/^g are slowly-varying profiles. As depicted in Fig. SI a and b, the momentum-space density profile of single realization 
of d opt = 500 and Q c = 2r shows the clear generation of these discrete ±nk p modes. Furthermore the Fourier transform of the 
coherence cr eg depicted in Fig. Sic shows clear superposition of ±nk p modes. 

The coherent SR-BEC interaction transfer most of the atoms to the ground state while the rest of them decay to other states 
due to the incoherent spontaneous processes as shown in Fig. Sid. Therefore the dynamics of the condensate after passing 
through the cavity is described by the single component Gross-Pitaevskii equation which reads 


l Ti 2 moj 2 z 2 - i 2 \ 

ihd t i// g - I H 2 f t* 88 1^1 j ^ g9 


(S 8 ) 
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FIG. SI. Panels a-d show the dynamics of the binary BEC inside the EIT cavity of d opt = 500 and Cl c = 2T within 2 /us. The evolution of 
the condensate densities in momentum space are shown in panels a and b where the high odd k p modes are generated in \jt g while the even k p 
modes are generated in xft e . Panel c shows the optical coherence cr eg . The coherent transfer of the BEC particle number from states \e) (green 
dashed-dotted line) to \g) (blue solid line) is shown. Panels e-h represent the time-of-fight simulation after the condensate is released from 
the cavity and harmonic trap. Panels e-g are the density components carry the momenta ±k p , ±3k p , and ±5k p , respectively. The total density 

profile \xjt g \ at t TO F = 1.8 ms. The contrast of the 2D plots is adjusted for better visualization. 


where we neglect the cross-species and atom-light interactions. 

During the time of flight (TOF) measurement, the condensate would split into several atomic clouds which corresponds to 
different ±nk p modes. To simulate the TOF dynamics, we numerically integrated equation. S8 by removing the trapping potential 
and interactions. As shown in Fig. SI e-h, we perform the TOF simulation for A t T 0 F = 1.8 ms. In Fig. SI e-g, the evolution of 
the condensed atoms carrying momenta are plotted and the velocities can be calculated from the slopes which agrees with the 
expected values, ±fik p /m , ±3 fik p /m, and ±hk p /m. In Fig. SI h, the total density distribution at t T 0 F = 1.8 ms is shown where 
six density bumps that are symmetric to the origin can be clearly observed. The two innermost density bump pair in Fig. SI h 
corresponds to the ±k p modes while the density bump pair located around z = ±0.5 mm is the ±3 k p modes and the outermost 
pair carries ±5k p momenta. 















